PROCEEDINGS 

OF SCIENCE 



o 
o 

O 

(N 

O 

Q 

On 



> 

m 

m 
(N 

o 



NLO QCD corrections to pp -> WW+jet+X 



Stefan Dittmaier*" Stefan Kallweif and Peter Uwer c t 

" Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut), D-80805 Milnchen, Germany 
h Faculty of Physics, University of Vienna, A- 1090 Vienna, Austria 

c Institutfiir Theoretische Teilchenphysik, Universitdt Karlsruhe, D-76128 Karlsruhe, Germany 



E-mail: Stefan . dittmaiergmppmu . mpg . de 


kallweit@mppmu . mpg . de 


uwer@particle . uni-karlsruhe . de 



We report on the calculation of the next-to-leading order QCD corrections to the production of 
W-boson pairs in association with a hard jet at hadron colliders, which is an important source 
of background for Higgs and new-physics searches. If a veto against the emission of a second 
hard jet is applied, the corrections stabilize the leading-order prediction for the cross section 
considerably. 
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1. Introduction 

The search for new-physics particles — including the Standard Model Higgs boson — will be 
the primary task in high-energy physics after the start of the LHC that is planned for 2008. The 
extremely complicated hadron collider environment does not only require sufficiently precise pre- 
dictions for new-physics signals, but also for many complicated background reactions that can- 
not entirely be measured from data. Among such background processes, several involve three, 
four, or even more particles in the final state, rendering the necessary next-to-leading-order (NLO) 
calculations in QCD very complicated. This problem lead to the creation of an "experimenters' 
wishlist for NLO calculations" ^ that are still missing for successful LHC analyses. The process 
pp — ► W + W~+jet+X made it to the top of this list. 

The process of WW+jet production is an important source for background to the production 
of a Higgs boson that subsequently decays into a W-boson pair, where additional jet activity might 
arise from the production or a hadronically decaying W boson [0]. WW+jet production delivers 
also potential background to new-physics searches, such as supersymmetric particles, because of 
leptons and missing transverse momentum from the W decays. Last but not least the process is 
interesting in its own right, since W-pair production processes enable a direct precise analysis of 
the non-abelian gauge-boson self-interactions, and a large fraction of W pairs will show up with 
additional jet activity at the LHC. 

In these proceedings we briefly report on our recent calculation [J|] of NLO QCD corrections 
to WW+jet production at the Tevatron and the LHC, but here we discuss results for the LHC only. 
Parallel to our work, another NLO study [Q] of pp — > W + W~+jet+X at the LHC appeared. A 
comparison of results of the two groups is in progress. 



2. Details of the NLO calculation 

At leading order (LO), hadronic WW+jet production receives contributions from the partonic 
processes qq — ► W + W~g, qg — ► W + W~q, and qg — > W + Wq, where q stands for up- or down- 
type quarks. Note that the amplitudes for q = u,d are not the same, even for vanishing light quark 
masses. All three channels are related by crossing symmetry. The LO diagrams for a specific 
partonic process are shown in Figure [I]. 

In order to prove the correctness of our results we have evaluated each ingredient twice us- 
ing independent calculations based — as far as possible — on different methods, yielding results in 
mutual agreement. 

2.1 Virtual corrections 

The virtual corrections modify the partonic processes that are already present at LO. At NLO 
these corrections are induced by self-energy, vertex, box (4-point), and pentagon (5-point) correc- 
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Figure 1: LO diagrams for the partonic process uu — > W + W~g. 
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Figure 2: Pentagon diagrams for the partonic process uti — * W + W~g. 

tions. For illustration the pentagon graphs, which are the most complicated diagrams, are shown in 
Figure || for a specific partonic channel. At one loop WW+jet production also serves as an off-shell 
continuation of the loop-induced process of Higgs+jet production with the Higgs boson decaying 
into a W-boson pair. In this subprocess the off-shell Higgs boson is coupled via a heavy-quark loop 
to two gluons. 

Version 1 of the virtual corrections is essentially obtained as for the related processes of ttH [j|] 
and tt+jet [§] production. Feynman diagrams and amplitudes are generated with FeynArts 1.0 
[^] and further processed with in-house Mathematica routines, which automatically create an out- 
put in Fortran. The IR (soft and collinear) singularities are treated in dimensional regularization 
and analytically separated from the finite remainder as described in Refs. [j|, ||]. The pentagon 
tensor integrals are directly reduced to box integrals following Ref. This method does not in- 
troduce inverse Gram determinants in this step, thereby avoiding numerical instabilities in regions 
where these determinants become small. Box and lower-point integrals are reduced a la Passarino- 



Veltman [10] to scalar integrals, which are either calculated analytically or using the results of 



Refs. [11]. Sufficient numerical stability is already achieved in this way, but further improvements 



with the methods of Ref. [12] are in progress. 



Version 2 of the evaluation of loop diagrams starts with the generation of diagrams and ampli- 
tudes via FeynArts 3.2 [|13|], which is independent of version 1.0 [{J. The amplitudes are further 



manipulated with FormCalc 5.2 [ 14] and eventually automatically translated into Fortran code. The 



whole reduction of tensor to scalar integrals is done with the help of the LoopTools library [M], 



which also employs the method of Ref. [|9p for the 5-point tensor integrals, Passarino-Veltman Q10| ] 



reduction for the lower-point tensors, and the FF package [15] for the evaluation of regular scalar 



integrals. The dimensionally regularized soft or collinear singular 3- and 4-point integrals had to 



be added to this library. To this end, the explicit results of Ref. []8j] for the vertex and of Ref. [ |16| ] 
for the box integrals (with appropriate analytical continuations) are taken. 

2.2 Real corrections 

The matrix elements for the real corrections are given by — » W + W~gggg and — ► 
W + W~qqq'q' with a large variety of flavour insertions for the light quarks q and q'. The par- 
tonic processes are obtained from these matrix elements by all possible crossings of quarks and 
gluons into the initial state. The evaluation of the real-emission amplitudes is performed in two 
independent ways. Both evaluations employ (independent implementations of) the dipole subtrac- 
tion formalism JT7| ] for the extraction of IR singularities and for their combination with the virtual 
corrections. 



Version 1 employs the Weyl-van-der-Waerden formalism (as described in Ref. [18]) for the cal 



culation of the helicity amplitudes. The phase-space integration is performed by a multi-channel 
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process 


a[pb] 


tfSherpa[pb] 


Aa/stat. error 


pp -> WW + ljet 


46.453(16) 


46.4399(94) 


+0.70 


pp -> WW + 2jets 


31.555(17) 


31.5747(63) 


-1.08 



Table 1: Comparison of LO cross sections with Sherpa (taken from Ref. [E4|) 



Monte Carlo integrator [19] with weight optimization [^] written in C++, which is constructed 
similar to Racoon WW [^l]]. The results for cross sections with two resolved hard jets have been 



checked against results obtained with Whizard 1.50 [22] and Sherpa 1.0.8 [23]. Details on this part 



of the calculation can be found in Ref. [24], the comparison to Sherpa results is briefly illustrated 
in Table [I]. 1 In order to improve the integration, additional channels are included for the integration 
of the difference of the real-emission matrix elements and the subtraction terms. 



Version 2 is based on scattering amplitudes calculated with Madgraph [ |25| ] generated code. The 
code has been modified to allow for a non-diagonal quark mixing matrix and the extraction of the 
required colour and spin structure. The latter enter the evaluation of the dipoles in the Catani- 
Seymour subtraction method. The evaluation of the individual dipoles was performed using a C++ 
library developed during the calculation of the NLO corrections for tt+jet [[J. For the phase-space 
integration a simple mapping has been used where the phase space is generated from a sequential 
splitting. 



3. Numerical results 



We consistently use the CTEQ6 Q26[ ] set of partem distribution functions (PDFs), i.e. we take 
CTEQ6L1 PDFs with a 1-loop running a s in LO and CTEQ6M PDFs with a 2-loop running a s 
in NLO. We do not include bottom quarks in the initial or final states, because the bottom PDF is 
suppressed w.r.t. to the others; outgoing bb pairs add little to the cross section and can be exper- 
imentally further excluded by anti-b-tagging. Quark mixing between the first two generations is 
introduced via a Cabibbo angle 6q = 0.227. In the strong coupling constant the number of active 
flavours is Nf = 5, and the respective QCD parameters are Aj° = 165 MeV and A^ s = 226 MeV. 
The top-quark loop in the gluon self-energy is subtracted at zero momentum. The running of a s 
is, thus, generated solely by the contributions of the light quark and gluon loops. The top-quark 
mass is m t = 174.3 GeV, the masses of all other quarks are neglected. The weak boson masses are 
M w = 80.425GeV, M z = 91.1876GeV, and M H = 150GeV. The weak mixing angle is set to its 
on-shell value, i.e. fixed by c 2 = 1 — s^, =M^/M^, and the electromagnetic coupling constant a 
is derived from Fermi's constant = 1.16637 x 10 _5 GeV 2 according to a = ^/IG^M^s^/k. 



We apply the jet algorithm of Ref. [27] with R = 1 for the definition of the tagged hard jet 



and restrict the transverse momentum of the hardest jet by /?T,jet > PTjet.cut- In contrast to the real 
corrections the LO prediction and the virtual corrections are not influenced by the jet algorithm. In 
our default setup, a possible second hard jet (originating from the real corrections) does not affect 



'The input parameters of Ref. [ p4| ] are set as below, apart from Os(Mx) = 0.1187 (1-loop evolved to the scale 
Mren = Mfact = Mw), and a CKM matrix in Wolfenstein parametrization (to 2nd order in A) with X = 0.22. The transverse 
momenta of additional jets are restricted by px.jet > 20GeV, and the jet-jet invariant mass by M(jet, jet) > 20GeV. No 
jet algorithm is applied, since genuine LO quantities are considered in Ref. |J24], 
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Figure 3: LO and NLO cross sections for WW+jet production at the LHC: scale dependence with renor- 
malization and factorization scales set equal to ji (lhs) and cut dependence (rhs) (taken from Ref. [B[]). 



the event selection, but alternatively we also consider mere WW+jet events with "no 2 nd separable 
jet" where only the first hard jet is allowed to pass the pjjet cut but not the second. 

The left-hand side of Figure || shows the scale dependence of the integrated LO and NLO 
cross sections at the LHC for pTjet,cut = lOOGeV. 2 The renormalization and factorization scales 
are identified here (/i = /i ren = ^f act ), and the variation ranges from }x = 0.1 A% to jX = 10 M\y. 
The dependence is rather large in LO, illustrating the well-known fact that the LO predictions can 
only provide a rough estimate. Varying the scales simultaneously by a factor of 4 (10) changes the 
LO cross section by about 35% (70%). 

Only a modest reduction of the scale dependence to 25% (60%) is observed in the transi- 
tion from LO to NLO if W pairs in association with two hard jets are taken into account. This 
large residual scale dependence in NLO, which is mainly due to gg-scattering channels, can be 
significantly suppressed upon applying the veto of having "no 2 nd separable jet". In this case the 
uncertainty is 10% (15%) if the scale is varied by a factor of 4 (10). The relevance of a jet veto in or- 



der to suppress the scale dependence at NLO was also realized [ ]28[ ] for genuine W-pair production 
at hadron colliders. 

Finally, we show the integrated LO and NLO cross sections as functions of pTjet,cut on the 
right-hand side of Figure |[ The widths of the bands, which correspond to scale variations within 
Mw/2 < /I < 2My/, reflect the behaviour discussed above for fixed value of /?T,jet,cut- 

For the LHC 

the reduction of the scale uncertainty is only mild unless WW+2jets events are vetoed. 
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